Disorder-induced Phase Transition of Vortex Matter in MgB 2 
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Measurements of single crystalline MgB2 with torque magnetometry in fields up to 90 kOe reveal 
a sharp peak in the irreversible torque at about 0.85 H C 2- In the region between peak onset and max- 
imum, pronounced history effects occur. Angle and temperature dependence of the characteristic 
peak fields track those of H C 2- The features observed suggest that the peak marks a disorder- induced 
phase transition of vortex matter between a quasi-ordered Bragg glass and a highly-disordered phase. 

PACS numbers: 74.25.Dw, 74.60.Ge, 74.70.Ad, 74.25.Ha 
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The new superconductor MgB2 is considered to have 
great potential for applications, and a lot oLresearch ac- 
tivity has concentrated on this compound.El Much light 
was shed on the superconducting mechanism, there is 
mounting evidence that MgB 2 is a two-band supercon- 
ductor with a substantial diffcrcnrjc-hctwcen the super- 
conducting gaps of the two bands.Emdlj About the super- 
conducting phase diagram, however, less is known. Most 
"phase diagrams" published contain only the upper criti- 
cal fields H C 2 . Here, the situation was clarified insofar as 
later bulk measurements all find a pronounced anii 
7 of H C 2, decreasing with increasing temperature,! 
although there are still discrepancies of the exact 7(T) 
dependences reported. MgB 2 is, particularly concern- 
ing the importance of thermal fluctuations and the value 
of k — A/£, intermediate between the high T c cuprates 
and low T c superconductors. Studying the vortex matter 
phase diagram of MgB 2 may thus help in understand- 
ing the phase diagrams of various superconductors in a 
unified way. 

From the study of cuprate superconductors is known 
that the H—T phase diagram contains more transition 
lines than the upper and lower critical fields. Identified 
were for example a melting transition between a quasi- 
ordered vortex lattice, called Bragg glass, and a disor- 
dered vortex fluidjilj as well as an order-disorder tran- 
sition between, the-Bragg glass and a hiahJy disordered, 
glassy phase O'HLil The latter transitiorO has been-ottj 
served also in low T c superconductors, such as Nb.Se 2 ,Ej't!3 
and even in the elementalpSuperconductor Nb,EZl but not 
in ultra-pure Nb crystalsJia This transition is generally 
associated with a peak in the critical current density j c 
and pronounced history effects. 

In single crystals of MgB 2 , a quasi-ordered vortex 
structure has been observed in low fields by scanning tun- 
neling spectroscopy,li3 showing that at least under some 
conditions a Bragg glass is the stable vortex phase. Since, 
by tuning the amount of quenched random point-like dis- 
order, the stabilization of a highly disordered phase can 
always be favored, an order-disorder transition in fields 
below H C 2 should be observable in MgB 2 as well, at least 
for certain impurity concentrations. Although a phase 
transition distinct from H c % or i? c2 has not been sug- 



gested yet in MgB 2 , a peak effect, and accompanying 
history effects have been observed in transport measure- 
ments for 7?||c.lI 

Here, we report the observation of a pronounced, sharp 
peak effect (PE) by torque magnetometry in fields close 
to, but clearly distinct from, H c %. A minor hysteresis 
loop (MHL) study shows pronounced history effects in 
the region between the onset and the maximum of the 
peak. Angle and temperature dependence of the charac- 
teristic fields is reported, and we propose a phase diagram 
for MgB 2 . 

The measurements were performed on a high-quality 
single crystal of MgB 2 , sample B of Ref. 0. The T depen- 
dence of the magnetization (upper inset of Fig. |l|) shows 
a sharp (0.3 K with a 10% — 90% criterion) transition 
to the superconducting state at 38.2 K, indicating a high 
quality of the crystal. Measurements to study the PE 
were carried out with the torque option of a Quantum 
Design 9T PPMS. Measurement runs consisted in vary- 
ing the field H at fixed angle 9 between H and the c-axis 
of the crystal, and recording the torque ~t — m x H , 
where m is the magnetic moment of the crystal. 

One of the curves measured is shown in the lower in- 
set of Fig. [l]. For better comparison with magnetization 
curves, t/H vs H is shown. The main panel shows a 
magnification of the PE region. The peak is well pro- 
nounced and very sharp. Various characteristic fields are 
indicated: The maximum of the peak for field increasing 
(if max) and decreasing (if max ) branch of the hysteresis 
loop, and the onsets of the peak, 7Jj n and H^ n . The sep- 
aration of the two onset fields is larger, similar to the 
case of the cuprate superconductors (see, e.g., Rcf. p0[ ). 
Also indicated is the irreversibility field i?i rr , where the 
two branches of the hysteresis loops meet. Thejaeak rc- 
sembles-qualitatively peaks observed in NbSe 2 ,E3'EJ and 
CeRu 2 .E3 

To investigate possible history dependences of j c , we 
performed several minor hysteresis loop (MHL) measure- 
ments in and around the peak: The field is cycled up and 
down by a small amount several times, ideally until the 
loops retrace each other, indicating that the vortex, sjis-. 
tern reached a stable pinned state in the given fieldJl3'E3 
MHL measured, within full loops, in four different regions 
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FIG. 1: Torque t/H vs field H at 14K and 77.5 deg. The 
direction of the field change is indicated by thick arrows. The 
irreversibility field H[ TI and the onset and maximum fields H on 
and f/max of the PE for the H increasing ( T ) and decreasing 
(■'•) branch are marked. Also shown are some of the MHL (see 
text) measured, labeled A-D. Upper inset: M(T) curve in the 
transition region, in a field H\\c of 1 Oe, zero field cooled (•) 
and field cooled (o). Lower inset: t/H vs H of the curve in 
the main panel, for the whole field range. 



of the PE are indicated in the figure (A-D). 

Torque t/H values of MHL A (Fig. 0a)) vary signifi- 
cantly as the MHL is cycled through repeatedly. Partly, 
this may be explained by relatively strong normal relax- 
ation processes. However, a pronounced difference can 
be seen between MHL started from the field increasing 
(H*) branch of the full hysteresis loop (FHL), and the 
one started from the field decreasing branch. The 

latter has a significantly higher width initially. This ef- 
fect can be explained by a difference in the vortex con- 
figuration between W and H* in the region of MHL 
A. In the configuration on H*, j c (proportional to the 
width of the MHL)c2lis higher, i.e. the vortices are pinned 
stronger. Repeated cycling causes the width of the MHL 
started from to approach the one started from W, 
indicating that the vortex configuration on is only 
meta-stable. History effects are even more pronounced 
for MHL B (Fig. ||b)). Here, the initial 7J T branch of 
the MHL started from the H l branch of the FHL (full 
line indicated by arrows) clearly is below the branch 
of the FHL (thick dashed), indicating larger hysteresis. 
This behavior contradicts Bean's critical state model,o 
where the hysteresis of partial hysteresis loops can never 
be higher than the one of the full loop. We point out that 
simple relaxation effects cannot account for this specific 
effect. It can be explained by the vortex configuration on 
the H* branch of the FHL (where the MHL was started) 
having a higher j c than the vortex configuration on the 
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FIG. 2: a),b),c) Magnification of MHL A, B, and D also dis- 
played in Fig. [y. MHL started from the field increasing branch 
of the full hysteresis loop are shown as dotted lines, while 
those started from the field decreasing branch are shown as 
full lines, d) Width of the hysteresis of MHL B started from 
the field increasing/decreasing (o/») branch of the full hys- 
teresis loop, as a function of cycling. 



H} branch. The variation of the hysteresis width with 
cycling (Fig. |2|d)) demonstrates the meta-stable nature of 
the vortex configuration on the H^- branch of the FHL, 
while the vortex configuration of the W branch of the 
FHL is stable, or close to. In contrast, no clear devia- 
tions in the cycling behavior between H* and branch 
started MHL are visible for MHL C and MHL D (Fig. 
§c)), as well as for a MHL measured in the region around 
68kOe (not shown). 

In summary, between H^ n and H^ ax , pronounced his- 
tory effects occur. They can be accounted for by the 
coexistence of a meta-stable high-field vortex configura- 
tion with high pinning and a stable low-field, low pinning 
configuration. Above -f/^ax and below H^ n , no significant 
history effects are observed, indicating that there is only 
one vortex configuration, which is stable. The larger hys- 
teresis width of MHL started from indicates pinning 
in the configuration stable above i?i ax to be stronger 
than pinning in the configuration stable below H* a . 

The variation of the peak onsets and maxima with ap-i 
gle at 18 K is shown, together with H c2 (0) and i2i rr (0),c3 
in Fig. [| Since the visibility of the peaks is diminished 
at higher temperatures, onsets and maxima were deter- 
mined from At(H) = t(H^) — t(H') curves. The char- 
acteristic peak fields follow the angular dependence of 
H C 2, as indicated by fits to the theoretical H C 2(9) de- 
pendence according to the anisotropic Ginzburg-Landau 
theory (see Ref. ||) , while the angular scaling of the irre- 
versibility field is less clear. This can be seen also in the 
inset, displaying the 9 dependence of the characteristic 
fields, reduced by the upper critical field. The onset field 
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FIG. 3: Angle dependence of various characteristic fields at 
18 K. Shown are the upper criticftli field H C 2 (0, from Ref. 
ph, the irreversibility field Hi TT (•),t3 the peak maximum field 
jHmax (a), and the peak onset field iinn (▼)■ Full lines are fits 
of the theoretical H C 2(8) dependences Dashed lines are guides 
for the eye. Inset: Angle dependence of reduced (divided by 
H C 2{9)) characteristic fields. 
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FIG. 4: Phase diagram of MgB2 single crystal at an angle 
of 77.5 deg between the c-axis of the crystal and the applied 
field: The temperature dependence of the characteristic fields 
H C 2, -Hmax and H on is given. They mark boundaries between 
the normal state and the various phases of vortex matter. The 
irreversibility field Hi TI is also shown. The inset shows the T 
dependence of the characteristic fields scaled by H C 2- 



is approximately constant at about 0.8 H C 2 and the max- 
imum field at about 0.85 H C 2- H m is located at about 
0.9i? C 2, but seems to get slightly lower as 9 — > 0. 

The characteristic fields could not be determined with 
enough accuracy in the whole region of angles: Since 
~t = m x H and m points, for H\\c or H\\ab, into the 
same direction, the sensitivity is much lower for angles 
close to and 90 deg. Due to the pronounced anisotropy 
of MgB 2 at 18 K (7 ~ 5.7)0 m tends to be directed al- 
most perpendicular to the planes, except at very high 
angles. Therefore, the maximum effective sensitivity of 
the torque magnetometer is achieved at angles in the re- 
gion of 75-80 deg. SQUID measurements performed on 
the same crystal with H\\c and iJ||a& showed no sign of 
a PE in the region around 0.8 H C 2- This is likely due to 
insufficient sensitivity of the SQUID and field inhomo- 
geneities in the SQUID magnetometer, which, due to the. 
movement of the sample, tend to smear sucb_features.c3 
However, preliminary ac susceptibility dataEH measured 
on the same crystal indicate the PE to be present both 
for H\\c and H\\ab, confirming that the underlying mech- 
anism is a feature for all directions of H . 

In Fig. the T dependence of the characteristic fields 
is shown, for 77.5 deg, which corresponds roughly to the 
angle where the PE is most visible. The peak amplitude 
is reduced quickly by increasing T, and above 22 K the 
PE is no longer clearly discernible in the t(H) data. This 
is due to the decreased sensitivity of the magnetometerEj 
and due to thermal smearing of the effective disorder po- 
tential. The inset shows that the positions of H on and 



therefore, unlikely that they would merge with the up- 
per critical field at some higher temperature. In contrast, 
H m shifts to lower fields relative to H C 2 as T increases, 
also likely due to a smearing of the effective pinning land- 
scape by thermal fluctuations. There is no indication 
that H{ rr may correspond to a phase transition. 

Before discussing the PE in terms of a disorder-induced 
phase transition of vortex matter, we briefly examine al- 
ternative origins of the PE. The possibility that the PE is 
due to inhomogeneities or extended defects is not likely. 
A second crystal grown with the same technique, but 
under slightly different conditions, shows also a PE, in 
similar fields. The pinning properties of the two crystals 
for H nearly aligned parallel, to the ab planes are dif- 
ferent in a pronounced wayHj It therefore seems rather 
unlikely that the two crystals would have identical struc- 
tural, non-intrinsic features leading to a similar PE. A 
further possibility would be a change of the elastic con- 
stants of the vortex lattice when iJr-approaches H C 2, not 
associated with a phase transition.^! However, the spe- 
cific form of the history effects observed in the PE region 
are hard to explain without a phase transition. Ther- 
mal melting can rather be excluded. Thermal fluctua- 
tions should be much less important in MgB2 than in the 
cuprate superconductors: MgB 2 has a Ginzburg num- 
ber Gi = |(7fc B T c / 'H^(0)^ b (0)) 2 , a measure of the im- 
portance of thermal fluctuations, of the order of_L0~ 5 
only, while the cuprates typically have Gi « 10 _2 .til On 
the other hand, thermal fluctuations should be morejm- 
portant than for example, in Nb with Gi w 10~ 10 ,E3 or 
NbSe 2 with Gi w 10~ 8 EJ A calculation of the melting 



ff max relative to H c i are approximately constant. It is, field H m , using Eq. (26) of Ref. &h leads, at 14 K, to 
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H m /H C 2 w 0.97, much higher than the location of the 
peak and therefore hardly can account for it J23 although 
it was showiL-that point disorder can shift H m to slightly 
lower fields.lEl Also, a liquid caused by thermal melting 
should have weaker pinning properties than the solid lat- 
tice. 

An important fact deduced from the MHL experiments 
is that the high field phase has got a higher critical cur- 
rent density than the low field phase. This is the case 
for the transition from a Bragg glass to a highly disor- 
dered glassy phase.Eil That we indeed observed this phase 
transition is supported by the pronounced history depen- 
dence of j c in the region between onset and maximum of 
the peak, of a form similar to observations of the PE in 
NbSe2 and not accountable for by relaxation effects. The 
location relatively close to H C 2 is expected for a supetj 
conductor with low Gi and relatively weak disorder .Eil 
In NbSe2, a superconductor with comparable, but even 
lower Gi, there is conclusive evidence that the PE is in- 
deed due to the transition between a Bragg glass and a 
highly disordered phase.Ej The history effects mark the 
region of meta-stability, where a macroscopic coexistence 
of the two phases is possible. Pinning of the phase bound- 
ary is directly responsible for the history effects. The lo- 
cation of the PE with respect to H C 2, together with the 
history effects studied and the observation of a higher 
critical current density in the high field vortex configura- 
tion, thus indicate that the PE in MgB2 marks the tran- 



sition between the Bragg glass and a highly disordered 
phase, which may be termed "amorphous" or "pinned 
liquid" . If the PE observed by Welp et alu is of the same 
origin, the larger separation of the PE from H C 2 in our 
case indicates that the crystal investigated by us has a 
higher amount of random point-like disorder. Further 
investigations of the transition line with controlled tun- 
ing of the amount of disorder, as was dope in the case 
of the cuprates, by electron irradiationt^l and chemical 
substitution,!^] may help finding a unified description of 
the phase diagrams of different superconductors. 

In summary, using torque magnetometry, we observed 
a pronounced, sharp peak effect in single crystalline 
MgB2. Onset and maximum of the peak are located at 
about 0.8 H C 2 and 0.85 H C 2, with little dependence on the 
temperature or the direction of the applied field. Peak 
form, history effects between onset and maximum, as well 
as the location of the peak are consistent with the peak 
effect marking a phase transition between the Bragg glass 
and a highly disordered phase of vortex matter. 
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Note After submission of this manuscript, we became 
aware of a reportE-3 on ac susceptibility measurements for 
H\\c, coming to similar conclusions. 



Email: angst@phys.ethz.ch 



For a review of early research see C. Buzea and T. Ya- 
mashita, Supercond. Sci. Technol. 14, R115 (2001). 

2 F. Bouquet et al, Europhys. Lett. 56, 856 (2001). 

3 H. Schmidt et al, Phys. Rev. Lett. 88 , 127002 (2002). 

4 M. Iavarone et al, |cond-mat/0203329 

5 M. Angst et al, Phys. Rev. Lett. 88, 167004 (2002). 

6 A. Sologubenko et al, Phys. Rev. B 65, 180505(R) (2002). 

7 S. L. Bud'ko and P. C. Canfield, Phys. Rev. B 65, 212501 
(2002). 

8 M. Zehetmayer et al, Phys. Rev. B 66, 052505 (2002). 

9 U. Welp et al, |cond-mat/0203337 . 

10 A. Schilling et al, Nature 382, 791 (1996). 

11 T. Giamarchi and P. Le Doussal, Phys. Rev. Lett. 72, 1530 
(1994); Phys. Rev. B 52, 1242 (1995); T. Nattermann and 
S. Scheidl, Adv. Phys. 49, 607 (2000). 

12 B. Khaykovich et al, Phys. Rev. Lett. 76, 2555 (1996); 
Phys. Rev. B 56, R517 (1997). 

13 D. Giller et al, Phys. Rev. Lett. 79, 2542 (1997). 

14 We follow convention and treat melting and order-disorder 
transition separately, although there is evidence that they 
form one line of phase transitions, which just changes 
character from thermally driven to disorder driven. See 
N. Avraham et al, Nature 411, 451 (2001). 

15 G. Ravikumar et al, Phys. Rev. B 63, 024505 (2001). 

16 M. Marchevsky, M. J. Higgins, and S. Bhattacharya, Na- 
ture 409, 591 (2001). 



17 
18 
19 
20 
21 
22 
23 

24 
28 



X. S. Ling et al, Phys. Rev. Lett. 86, 712 (2001). 

E. M. Forgan et al, Phys. Rev. Lett. 88, 167003 (2002). 

M. R. Eskildsen, private communication. 

M. Angst et al, Phys. Rev. B 65, 094518 (2002). 

V. V. Eremenko et al, Low Temp. Phys. 28, 6 (2002). 

S. B. Roy et al, Phys. Rev. B 62, 9191 (2000). 

The form of the MHL clearly indicates that a full reversal 

of the critical stateEj is achieved. 

C. P. Bean, Phys. Rev. Lett. 8, 250 (1962). 

A criterion T ltT /H sin(2#) = 10 _6 dyncmOe _1 was used for 

defining i?j rr . 

G. Ravikumar et al, Physica C 276, 9 (1997). 
R. Puzniak et al, in preparation. 

The signal-to-noise ratio is decreased in lower fields, since 
r oc H. 

M. Angst et al, |cond-mat/0206407 . 

A. B. Pippard, Philos. Mag. 19, 217 (1969). 

G. P. Mikitik and E. H. Brandt, Phys. Rev. B 64, 184514 

(2001). 

Combined ac susceptibility and dc magnetization measure- 
ments found the PE maximum to be located in fields equal 
to or even higher than the magnetization step signifying 
melting. See e.g. J. Shi et al, Phys. Rev. B 60, R12593 
(1999); T. Ishida et al, ibid 56, 5128 (1997). 
T. Nishizaki et al, Phys. Rev. B 61, 3649 (2000). 
M. Pissas et al, Phys. Rev. Lett. 89, 097002 (2002). 



